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Abstract 



Statistical features of "bursty" behaviour in charged and neutral fluid 
turbulence, are compared to statistics of intermittent events in a GOY shell 
model, and avalanches in different models of Self Organized Criticality (SOC). 
It is found that inter-burst times show a power law distribution for turbulent 
samples and for the shell model, a property which is shared only in a par- 
ticular case of the running sandpile model. The breakdown of self-similarity 
generated by isolated events observed in the turbulent samples, is well repro- 
duced by the shell model, while it is absent in all SOC models considered. On 
this base, we conclude that SOC models are not adequate to mimic fluid tur- 
bulence, while the GOY shell model constitutes a better candidate to describe 
the gross features of turbulence. 
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PACS. 05.45.-a: Nonlinear dynamics and nonlinear dynamical systems 
PACS. 05.65.+b: Self-organized systems 
PACS. 52.35.Ra: Plasma turbulence 
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Turbulence in neutral and charged fluids exhibits bursty behaviour, power laws and 
1// spectra. Understanding the origin of the rich and complex dynamics underlying these 
properties, is a challenging and fascinating task. In the years different models have been 
proposed aimed to catch the fundamental mechanism of this phenomenology. Among them 
simplifled models (" sandpiles" ) based on Self-Organized Criticality (SOC) |I| and dynamical 
systems approach 0, are particularly suitable for this purpose. SOC has been proposed as 
a paradigm for some complex systems, such as earthquakes, solar flares, river floods, fusion 
plasmas, etc. Despite the very simple governing rules, these systems exhibit a rich dynamics 
showing bursty behaviour with impulsive events ("avalanches" in sandpiles) which remind 
a variety of natural phenomena, including turbulence. Power laws are observed for the 
total energy |l|,|5J^ , the duration and the peak energy distributions of the bursts and for 
the total energy frequency spectra On the other hand models based on the dynamical 
systems approach have been developed by reducing the fluid equations, in order to isolate 
the basic mechanisms generating turbulence. An example of a model belonging to this class 
is the so called GOY shell model (see |0 and ref.s therein). This model describes a one 
dimensional fluid system by a set of ordinary differential equations of complex variables non 
linearly interacting through the coupling of the neighbour and next-neighbour modes. 

To what extent these simplifled models are good candidates to mimic turbulence is an 
intriguing question, reported, for example, in Ref.s |jl|J§-H]. The present paper is aimed to 
contribute to clarify this issue by comparing the predictions of these different dynamical 
models with samples of turbulence in neutral and charged fluids. The comparison is carried 
out by discussing two basic properties of turbulence, namely the presence of long time 
correlations and intermittency. For this purpose two statistical analysis tools have been 
used: the Distribution of the Waiting Times (WTD) At between subsequent bursts ||^ and 
the comparison of the Probability Distribution Function (pdf) of fluctuations resolved at 
different time scales p. These methods have been applied to SOC and GOY models as well 
as to experimental data. In particular, concerning SOC models, we have analyzed succesive 
generalization of the original sandpile model (BTW model), made by Kadanoff et al. 
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which varied the underlying microscopic rules and by Hwa and Kardar [|T0| which included 
finite drive effects (the running sandpile). 

It is worth noticing that the analysis of dynamical properties of these models requires a 
correct definition of the timescale, a point often neglected in SOC literature. For example 
Ref. 10 does not give any definition of the timescale, since it is focused on integral properties 
of avalanches, namely the total released energy and particles. With this in mind we choose 
to introduce in Kadanoff sandpile the same timescale definition proposed in Ref. . Accord- 
ingly in this case we built the signal of dissipated energy exit) considering every avalanche 
as an istantaneous event, with an amplitude equal to the integral avalanche size and the 
time between two subsequent events proportional to the number of particles needed to put 
the system in the unstable state again. On the other hand in the case of the running sanpile 
T0| the dissipated energy signal and its timescale are clearly defined. In each temporal step 



the system is continuously fed at random with a certain finite deposition rate probability 
Jin and the unstable sites are simultaneously updated. The dissipated energy signal enit) 
in this case is taken as the total number of unstable sites at each step. The timescale of 
the previous model is recovered in the limit Jj„ — 0. In our simulation we consider a sand- 
pile of size i, and associate to each lattice site n the height of the pile h{n,t). If there is 
any unstable site according to the rule h{n,t) — h{n ± > A, the system is evolved by 
taking a fixed amount Nf of sand from the unstable n site and redistributing grains to the 
neighbouring nil site (we set the threshold A = 50 and Nf = 20). 

The WTD -P(At) for the Kadanoff model (right-hand panel of Fig. |^) shows a clear 
exponential decay. This behaviour, which is known to be associated to independent Pois- 
sonian events 0, indicates that no long term correlations are present between subsequent 
avalanches, as found in the standard BTW model Q]. 

In the case of running sanpile at small J^, an exponential WTD was found (top right- 
hand panel of Fig. |^), recovering the behaviour of the Kadanoff sanpile. On the other hand, 
increasing the drive Jj„ a tendency to a power law behaviour P{At) ~ At"^-^ emerges. In 
order to measure the waiting times, we note that while at Jm avalanches are isolated 



events so that waiting times coincide with quiescent intervals, at high Jj„ there is a continuos 
activity. Thus to define the waiting times between subsequent events a threshold has been 
chosen applying the algorithm reported in [Q. The waiting time corresponds to the time 
interval between two subsequent events with amplitude beyond this threshold. In the limit 
Jin this method correctly gives threshold 0. 

It is known that long range time correlations can be obtained by overlapping signals from 
uncorrelated avalanches In order to test the origin of the power law above reported, 

following a procedure originally proposed in ref. a signal has been generated by randomly 
overlapping isolated avalanches for the case Jj„ —>■ 0, keeping an average dissipation rate 
equal to that for the case Jj„ = 4. In this case an exponential decay of the WTD has been 
found. Thus we conclude that the origin of power law in the WTD of the Hwa-Kardar 
model is due to the interaction among the avalanches inside the system. The emerging of a 
power law behavior with exponent close to /? ~ 1.5 in the WTD is a common feature found 
for all Hwa-Kardar systems, indipendently from their size i. 

In order to characterize the properties of the dynamics of these SOC systems, the energy 
fiuctuations 6er = e(t + r) — e(t) at different scales r of the energy signal e(t) defined above 
have been analysed. The increments (5e^ of a stochastic process are self-similar with respect 
to the scaling transformation r — Ar (with A > 0) if 6er ~ X^'^Sexr 0, which means that the 
cumulative distributions of the first and the second term coincide for a given scaling exponent 
a. This can be translated on the relationship between density probability of increments at 
two different scales r, namely ]9c//((5e^) = X^°'pdf{X^°6exr ). Test for self-similarity of random 
processes can be made without actually assuming a value for a. In fact it can be shown that 
pdfs of rescaled fluctuations 5Er = 5er/(J^{T)^ where (t^{t) = {([e(t + T) — e(i(:)]^)}^/^ (brackets 
being time averages), are invariant in a self-similar process, say pdf{6ET-) = pdf^SExr) 0- 

The pdfs of the fluctuations 6Et- at different scales r for the Kadanoff sandpile are shown 
on the left-hand panel of Fig. |1|. The bursty appereance of SOC signals is reflected onto 
the shape with fat tails of the pdfs. However pdfs are invariant when plotted in the rescaled 
form, so the process turns out to be purely self-similar. 
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Concerning the running sandpile, in the case of small drive {Jin — > 0) the result turns 
out to be close (upper-left panel of Fig. to the Kadanoff model, confirming that in the 
above limit the two models are equivalent. 

In the case of finite drive (bottom-left panel of Fig. the shape of pdfs approaches 
a Gaussian shape due to the interaction of the avalanches, and, even more noticeable, the 
underlying process remains purely self-similar. These results are consistent with the intrinsic 
self-similar structure of SOC models [ll], which is at the origin of its power law distributions 
and scaling properties It is worth noticing that self similarity is observed even in the 
case of finite drive, despite the common belief that " the true SOC behaviour arises only in 
the limit that the forcing rate approaches 0" [l^jlll . 



The MHD COY shell model describes the dynamic of the energy cascade in tur- 



bulence. The wavevector space is divided in shells, each shell being characterised by a 
discrete wavevector kn = 2"/co {n = 0,1, N), and by the (complex) dynamical variables, 
velocity M„(t) and magnetic field bn{t), which represent characteristic fluctuations across 
eddies at the scale ~ k~^. The dynamical behaviour of the model is described by the 
following set of ordinary differential equations 



dt 

1 1 -| * 

-j{Un-lUn+l-hn-lhn+l)--{Un-2Un-l-hn-2hn-l)] +fn (l) 



lyk^Un + ikn[{ 



^ = -vklbn + ^|(m„,+i6„+2 - bn+iUn+2) 

dt 
+ {Un-lbn+l - bn-lUn+l) + (u„_26n-l " &n-2^n-l)} (2) 

u and 1] are the kinematic viscosity and the resistivity respectively, and /„ is an external 
forcing term. The nonlinear terms have been obtained by imposing quadratic nonlinear 
coupling between neighbouring shells and the conservation of three ideal invariants. Detailed 
informations on the model and on simulations can be found in Ref.s P,p!2[|. 



The model has a fixed point in the form of a Kolmogorov-like scaling m„ ~ 6„ ~ k^^^^, 
and exibits a chaotic dynamics on a strange attractor in the phase space around the fixed 

6 



point. This is due to the presence of the ideal invariants. What is interesting here is the fact 
that fluctuations of the degree of chaoticity in the model generate a multifractal probability 
measure on the attractor. As a consequence the statistics of a sample of turbulence coming 
from simulations of the system is characterised by multifractal scalings, and this manifest 
itself as a breakdown of self-similarity of the system. That is, at variance with sandpile 
models, pdfs of the rescaled variables 5un = Real{un} / <Ju„{^n) and 6bn = Real{bn} / crtni^n) 
(see figure ||) strongly depend on the scale £„. At large scales pdfs are gaussian, while they 
tend to develop fat tails for small £„, towards the dissipative range. The departure from 
gaussianity visible at smaller scales, is due to bursts localised in time, which can be identified 
and extracted through an iterative procedure on and |6nP; described extensively in Ref. 
0. The WTD for these events is shown in the bottom panels of figure |^. We found a power 
law -P(At) ~ At~^ over at least two decades (the scaling exponents /? are reported on the 
figures). Long-range correlations between bursts are the origin of the power law for waiting 
times. 

We proceed now to investigate some samples of turbulent flows. First of all we con- 
sider a sample of turbulence as measured in the atmospheric boundary layer [|r^, where 
the longitudinal velocity field u{t) is recorded. Then we consider two samples of magnetic 
field intensity B{t) in plasmas as recorded both in a given low-speed stream of the solar 
wind by Helios II spacecraft [|14|, and in a laboratory device for thermonuclear fusion re- 



search (the RFX experiment |T^). Even if turbulent samples are collected in very different 
physical conditions, they show the same statistical properties. The rescaled fluctuations, 
namely 6ur = [u{t + r) — ■u(t)]/(T,u(r) and 6Br = [B{t + r) — B(t)]/(TB{T), depend on the 
scale T (figure ^). At large scales pdfs are gaussian, while they tend to develop fat tails as 
r — » 0, towards the dissipative range. The departure from gaussianity, which results in a 
correction to the usual Kolmogorov scaling laws, is due to " intermittency in fully developed 
turbulence" (see Ref. PJT6[1 and references therein). Since at smaller scales the wings of pdfs 
increase, intermittency is represented by intense and localised fluctuations which accumu- 
lates at smaller scales. These "events" are generated during the turbulent energy cascade. 



and can be interpreted as due to a non homogeneous breakdown of " eddies" towards small 
scales Their presence originates the departure from a pure self-similarity of fluctua- 
tions at the various scales. The time of occurrence of the events in a turbulent flow can be 
extracted through the procedure described in Ref . , and the distribution of waiting times 
between them is shown in the bottom panels of flgure ^ for the different samples. A power 
law over two decades, with non universal scaling exponents /?, is found in all samples. 

To summarize, we found that turbulence and the GOY MHD shell model always show 
power law for WTD, while this property is shared only by a SOC model, namely the running 
sandpile in the "strong" feeding regime, far from BTW limit p],p]T| . Moreover pdfs of rescaled 
fluctuations in both BTW and other SOC models do not change shape with the scale r. On 
the contrary intermittency in turbulent flows manifests itself as a breakdown of a pure self- 
similarity, leading to increasingly non gaussian tails for pdfs as the scale becomes smaller. 
The same phenomenon is observed in the GOY shell model, and is due to "time intermit- 
tency" . Non linear terms in the model, which conserve some few invariants, are responsible 
for the occurrence of bursts of chaoticity which are concentrated on the dissipative shells 
Such typical behaviour in the shell model generates a breaking of the global scaling 
invariance of the model, and this results in a correction to the Kolmogorov-like scaling law. 
This correction has the same statistical effects of the corrections of the Kolmogorov scaling 
laws due to usual intermittency in real fluid flows. Then as far as the statistical properties 
we investigated are concerned, SOC models we examined do not adequately mimic neither 
fluid nor plasma turbulence. On the contrary time intermittency in the GOY shell model 
allows the model to capture the gross statistical features of real turbulence. 
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FIG. 1. We show the results obtained from the KadanofF "local limited" sandpile obtained with 
a system size £ = 1024. In the left panel we show the pdfs of the rescaled fluctuations at different 
scales r, and in the right panel (log-lin plot) we show the distribution of waiting times between 
avalanches. Simulations obtained with different values of £ show the same results. 
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FIG. 2. Results obtained from the running sandpile for two different values of Jin (see text). In 

the left panels we show the pdfs of the resealed fluctuations at different scales r, and in the right 
panels (log-log plot) we show the distribution of waiting times between avalanches. 
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FIG. 3. We show results obtained from simulations of the GOY MHD shell model. Left-hand 
and right-hand columns refer to velocity and magnetic variables respectively. In the first three 
panels of each column we show the pdfs of the rescaled dynamical variables at three different scales. 
In the bottom panels of each column we report the distribution of waiting times between events at 
the scale £n/L = 2~^° (where L = k^^). The values of P represent the best fits of the power laws. 
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FIG. 4. Results obtained from a sample of fluid turbulence (left-hand column), a sample of 
magnetic turbulence measured by the Helios II satellite in the interplanetary space at distance 
R = 0.9 Astronomical Units (central column), and a sample of magnetic turbulence measured in a 
magnetically confined plasma in Reversed Field Pinch configuration (right-hand column). In the 
first three panels of each column we show the pdfs of the rescaled dynamical variables at three 
different scales. In the bottom panels of each column we report the distribution of waiting times 
between events at the smallest scale. The values of (3 represent the best fits of the power laws. 
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